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Itinerant Ferromagnetism and p + ip' Superconductivity in Doped Bilayer Silicene 
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School of Physics, Beijing Institute of Technology, Beijing 100081, China 

We study the electronic instabilities of doped bilayer silicene using the random phase approxima¬ 
tion. In contrast to the singlet d + id' superconductivity at the low doping region, we find that the 
system is an itinerant ferromagnet in the narrow doping regions around the Van Hove singularities, 
and a triplet p + ip' superconductor in the vicinity of these regions. Adding a weak Kane-Mele 
spin-orbit coupling to the system further singles out the time-reversal invariant equal-spin helical 
p + ip' pairing as the leading instability. The triplet pairing identified here is driven by the ferro¬ 
magnetic fluctuations, which become strong and enhance the superconducting critical temperature 
remarkably near the phase boundaries between ferromagnetism and superconductivity. 
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I. INTRODUCTION 

Magnetism and unconventional superconductivity 
(SC) as well as the intimate interplay between them 
have been the focuses of condensed matter physics for 
decades due to their rich physics and important applica¬ 
tions. Among these subjects, the realizations of itinerant 
ferromagnetism (FM) and triplet SC are of particular 
importance in recent years. In general, the triplet SC 
[I], which is connected with topological SC [2, 3 and 
becomes hot topic recently, is believed to be driven by 
ferromagnetic spin fluctuations near the FM order. How¬ 
ever, the realization of itinerant FM from the Stoner cri¬ 
terion |4 usually requires finite and, most of the time, 
strong electron interactions [SHIS], which is hard to deal 
with in the weak coupling perturbative approaches. One 
way to overcome this difficulty is through the introduc¬ 
ing of divergent density of states (DOS) at the Van Hove 
(VH) singularities of the system, which can induce these 
instabilities without strong electronic interactions. It’s 
proposed recently that, for a system with its Fermi sur¬ 
face (FS) doped to time-reversal (TR) variant VH saddle 
points, weak repulsive electron interactions can usually 
drive itinerant FM and triplet SC M- 

On another front, as the Si-based counterpart of 
graphene, silicene has been synthesized recently mm, 
with experimental evidence showing possible SC in the 
doped case [20] . Furthermore, bilayer silicene (BLS) has 
also been available 21 . with the energetically most fa¬ 
vored stacking way between its two layers identified by 
first-principles calculations [22]. Based on the metallic 
band structure of undoped BLS, the antiferromagnetism 
(AFM) and the chiral d j- id' SC tuned by strain have 
been proposed [22] . This intriguing result motivates us 
to further investigate the electronic instabilities in doped 
BLS, specifically focusing on the VH doping levels since 
the divergent DOS there favors the occurrence of elec¬ 
tronic instabilities. Paying our attention to VH doping, 
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we notice that, in VH-doped monolayer graphene whose 
VH saddle points locate at TR invariant momenta, the 
chiral spin density wave or the chiral d + id' pairing has 
been proposed [231427] . Similar results have also been 
found in monolayer silicene [28] . In contrast, the inter¬ 
esting property of the VH singularities here in BLS lies in 
that the VH saddle points locate at TR variant momenta. 
For such VH singularities, the study based on the renor¬ 
malization group theory has pointed out the possibility 
of the formation of itinerant FM and triplet SC El- 

In this paper, we perform the calculations based on 
the random phase approximation (RPA) to investigate 
possible electronic instabilities of doped BLS. The main 
results of our calculations are as follows. In addition to 
the d + id' SC occurring at low doping levels, the itiner¬ 
ant FM and the triplet p + ip' SC emerge as the leading 
instabilities of the system in the narrow doping regions 
around the VH singularities and the vicinity of these re¬ 
gions respectively. In the presence of a weak Kane-Mele 
spin-orbit coupling (SOC), the helical p + ip' pairing be¬ 
comes the leading instability of the system. The emer¬ 
gence of the FM and the triplet SC results from the large 
DOS and the strong ferromagnetic correlation near the 
VH singularities. Near the critical doping level separat¬ 
ing the FM and triplet SC, the strong ferromagnetic fluc¬ 
tuations will greatly enhance the superconducting criti¬ 
cal temperature, which provides possibility to realize this 
triplet p + ip' pairing state at experimentally accessible 
temperatures. 

The rest of this paper is organized as follows. In Sec. 
II, we describe the Hubbard model of BLS, as well as 
the RPA approach. In Sec. Ill, we calculate the sus¬ 
ceptibilities of the system, and demonstrate the itinerant 
FM occurring around the VH singularities. In Sec. IV, 
we study the superconducting pairing symmetries for dif¬ 
ferent doping levels, and propose that the p + ip' pairing 
dominates over the d + id' one in the vicinity of the ferro¬ 
magnetic regions. Finally in Sec. V, a conclusion will be 
reached after discussions on the experimental detection 
of the novel p + ip' pairing state proposed here. 
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II. MODEL AND APPROACH 

The lattice structure of BLS [22] is shown in Fig. [lja), 
which belongs to the D^d point group. While sublattice 
A\ of the upper silicene layer couples vertically to sub¬ 
lattice A 2 of the lower layer with a bond-length l v = 2.52 
A, the two sublattices Ai and Bi within the same layer 
l (= 1,2) couples to each other with a bond-length 
l n = 2.32 A. Approximately equal bond lengths, together 
with the bond-angle 0 = 106.65° between the two bonds 
describes an orbital hybridization more like the sp 3 type 
than the planar sp 2 one. This lattice structure leads 
to a strong interlayer coupling, and the resulting strong 
bonding-antibonding splitting between orbitals A 1 and 
A 2 pushes them far away from the Fermi level. Thus 
the low energy subspace formed by orbitals B\ and B 2 
will take responsibility for the main physics of the system 
[22] . This feature of BLS is obviously different from that 
of bilayer graphene. 



FIG. 1. (a) Optimized lattice structure of BLS. (b) The cor¬ 
responding band structure. In (a), both the top view (upper) 
and side view (lower) are shown. The intralayer nearest neigh¬ 
bor bond length Z n , the vertical bond length l v , and the angle 
6 between them are marked, together with the hopping inte¬ 
grals t n , ti, £ 2 , and £ 3 . In (b), the black dashed, red and blue 
solid horizontal lines denote the Fermi levels of the undoped, 
electron and hole VH doped systems respectively. 

According to Ref. [22] . the low energy physics of BLS 
near the FS can be described by the following 4-band 
Hubbard model of the system: 

H ^ ^ C'j^ ao .H a /d(^k^)Cp i /d cr T U ^ ^ (1) 

kcrotf3 ia 

Here cr, a (/?), and i denote the spin, orbital, and unit cell 
indices respectively, and H(k) is the 4-band tight-binding 
(TB) Hamiltonian in the basis {|F>i}, \B 2 ),\Ai),\A 2 )}. 
The explicit expression of the TB Hamiltonian reads [22] 
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Here f(k) 

- e lk Ra with 

R a (a = 1,2,3) being 

the nearest-neighbor 

vector, A 

= —0.069 eV is the ef- 


fective on-site energy difference between atoms A[ and 
H/, the hopping integrals £ n = 1.130 eV, t\ = 2.025 
eV, t 2 = 0.152 eV, and £3 = 0.616 eV. Since the basis 
(|F?i), \B 2 ), \A\), \A 2 )} is mainly composed of the 3 p z or¬ 
bital of silicon [29] , we set U = 1 eV as a rough estimate 
of the Hubbard interaction. 

The band structure for the above TB Hamiltonian ([2| 
is shown in Fig. 0 b). One feature of the band struc¬ 
ture is the 300 meV overlap between the valence and 
conduction bands near the K-points. For the undoped 
case, this overlap causes six pairs of small electron and 
hole pockets around and near the K-points respectively as 
shown in Fig. (2|a). The undoped system is thus intrinsi¬ 
cally metallic and can enter a superconducting state [ 22 ]. 
When the system is doped, regardless of electron or hole 
doping case, the shape of the FS grows gradually from 
separated electron and hole pockets first to six big merged 
pockets around the K-points ^b) and[2jd)), which fi¬ 
nally connect to one another at the VH saddle points, 
causing the Lifshits transition of the FS Qc) and 0 e)). 
Defining the doping level by x = n e — 1 where n e is the 
number of electrons per site, we find the doping level for 
the VH singularities are x = 0.2345 for electron doping 
and x = —0.1861 for hole doping respectively. The Fermi 
level for these VH dopings are marked in Fig. [TJb) , where 
the flatness of the band structure on the FS leads to the 
logarithmically divergent DOS near the Fermi level, as 
shown in Fig. We shall focus on these VH dopings 
in the following study because the divergent DOS around 
there urges the formation of itinerant FM and the result¬ 
ing strong ferromagnetic fluctuations in the vicinity of 
the FM region will induce high-temperature triplet SC. 

A special feature of the VH singularities of this mate¬ 
rial lies in that its VH saddle points locate on the M-T 
axes rather than at the TR invariant M-points as in bi¬ 
layer graphene as well as monolayer graphene or silicene. 
Such TR variant VH saddle points are named as “type- 
II” VH saddle points in Ref. [14], in contrast to the TR 
invariant “type-I” VH saddle points. The “type-II” VH 
singularities is special in that it allows for the formation 
of triplet SC. If the FS of a system contains TR invari¬ 
ant “type-I” VH saddle points, the triplet pairing will 
not be energetically favored because its odd parity gap 
function will have nodes at these TR invariant VH mo¬ 
menta, which is no good for the energy gain. On the 
contrary, the TR variant “type-II” VH saddle points of 
BLS locating on the M -Y axes provide the possibility for 
the system to enter the triplet pairing state. 

To study the electron instabilities of this system repre¬ 
sented by the Hubbard model (Hi), we adopt the standard 
multi-orbital RPA approach [ 22 J f30Tf34] . We first define 
and calculate the bare susceptibility tensor (< 7 , t). 

After that, the renormalized spin(s) or charge(c) suscep¬ 
tibilities xtu hM (*L r ) are °bt a i ne d in the RPA level. 
For each doping level, there will be a critical interaction 
strength Z7 C . For repulsive U > U c , the renormalized spin 
susceptibility diverges, implying the formation of long- 
range magnetic order. For U <U C , through exchanging 
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the spin or charge fluctuations, we obtain the effective 
pairing potential V af3 (k,q). Solving the linearized gap 
equation for V a P(k, q ) as an eigenvalue problem, we ob¬ 
tain the leading pairing gap function as the eigenvector 
corresponding to the largest eigenvalue. 






FIG. 2. (a)-(e) The shapes of the FS for different doping 

levels. The black dots in (c) and (e) indicate the VH saddle 
points, (f) The doping dependence of the DOS near the FS. 
The vertical dashed lines indicate the VH singularities. 


III. ITINERANT FERROMAGNETISM 


The bare (U = 0) susceptibility tensor of the model is 
defined as 

X { st Pq {k,T) = pY <r T 4(fci,r)cg(fei + fe,r) 

kik 2 

c\{k 2 + fc,0)c t (fc 2 ,0)) o , (3) 

Here (• • • )o denotes the thermal average for U = 0, T r 
denotes the time-ordered product, and p, g, s, t = 1, • • • ,4 
are the sublattice indices. Fourier transformed to the 
imaginary frequency space, the bare susceptibility can 


be expressed by the following explicit formulism, 


x ^ pq {k,iu n ) y, + fc ) 


k'a(3 


£?*(■ k’ + k)- 




icon + - ei, 


( 4 ) 


" k'+k 


Here iuj n is the Matsubara frequency, cq (3 = 1, • • * ,4 are 
the band indices, rip is the Fermi distribution function, 
^ and i a {k) are the eigenvalue and eigenvector of the 
TB Hamiltonian (J2|. The Hermitian static susceptibility 

matrix is defined as Xp°s(&) = iuj = 0). For each 

fc, the largest eigenvalue of this matrix represents 

the static susceptibility of the system in the strongest 
channel, and the corresponding eigenvector describes the 
pattern of the dominant intrinsic spin correlation in a 
unit cell of the system. 
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FIG. 3. (a)-(e) The k- space distributions of the zero- 

temperature static susceptibility x^(fe) for different doping 
levels. Typical (f) ferromagnetic pattern at the VH doping 
and (g) antiferromagnetic pattern at zero doping in a unit 
cell of the system. 


In Figs. Ia)B e), we show the fc-space distributions 
of the zero-temperature static susceptibility for 

different doping levels, which reveal the doping evolution 
of the static susceptibility. In particular, when the dop¬ 
ing level changes gradually from zero to the VH doping, 
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regardless of electron or hole doping case, the momenta 
of the maximum susceptibility evolve from the T-point 
(|3ja)) first to the points around it p|b) and [3^d)), and 
finally back to the T-point again (pFc) andpfe)). Such 
a doping evolution of the susceptibility originates from 
the evolution of the FS mentioned before, and indicates 
that the intra-sublattice spin correlation of the system 
changes gradually with doping from ferromagnetic first 
to antiferromagnetic, and finally back to ferromagnetic 
again. 

From the eigenvector corresponding to the largest 
eigenvalue of Xp°s(fe), we finds that the spin correlation 
within a unit cell is ferromagnetic-like (see Fig. [3jf)) near 
the VH doping levels, and antiferromagnetic-like (see Fig. 
[3^g)) near zero doping. Therefore, although the intra¬ 
sublattice spin correlations in both the VH doped and 
undoped systems are ferromagnetic, the inter-sublattice 
spin correlations in the former and latter cases are ferro¬ 
magnetic and antiferromagnetic respectively. 

When the interaction is turned on, we define the charge 
(c) and spin (s) susceptibilities for the model as 

xit )Pg (k,r) = T Y ( T T c bi( fc i’ r ) 

k 1 k 2 a 1 a 2 

Cqa, (fcl + k, r)c\ a2 (fc 2 + k, 0 )c t<T2 (fc 2 ,0)> , (5) 

xil )pq (k,T) = T Y 0-1(72 (T T 4 CTi (fei,r) 

k\k 2 cr \<j 2 

Cq ai (ki + k,r)c\ a2 {k 2 + k,0)c trT2 (fe 2 ,0)> , (6) 
xit )Pq {k,r) = -t Y (Trcl t (ki,r) 

k\ k 2 

Cqi{ki + k, r)4 i (fc 2 + fe, 0)c n (fc 2 ,0)) , (7) 

xit +)Pq (k,r) = F Y ( T A^i,r) 

k\ k 2 

c 9t (fei + fe, r)4 t (fc 2 + fe, 0)c u (fc 2 ,0)) , (8) 


where <ti, g^ =t? i are spin indices. For nonmagnetic 
states, we have x^ = x^ s+ ^ = X^ s = X^- For 
U — 0 , we further have x ^ = X ^ = X®- 

In the standard RPA approach [ 22 ; 301134] , the charge 
(spin) susceptibility for the model is given by 


X {c(s)) (k,ioj n )= I ± x ( 0 ) (fc, iio n )(U) x i0 \k,ioj n ) 


-1 


( 9 ) 


where (U) is a 16 x 16 matrix, whose only four nonzero 
elements are (U)w = U (fi = 1, • • * ,4) [ 22 ]. Clearly, 
the repulsive Hubbard interaction here suppresses x^ 
and enhances x^. When the interaction parameter U is 
weak enough, the RPA works well since all eigenvalues of 
the denominator matrix [.I db x® (fe, iuj n )(U)\ in Eq. ^ 
are positive and hence the matrix itself has an inverse. 
However, if U exceeds a critical value U c at which the 
lowest eigenvalue of [I — x^(fe,)(E/)] touches zero, 


the renormalized spin susceptibility %^ would diverge, 
which implies the formation of long-range magnetic or¬ 
der. 

The doping dependence of the critical interaction 
strength U c is shown in Fig. [4] The most obvious fea¬ 
ture of Fig. [4] is that U c drops abruptly to zero near the 
VH singularities due to the divergent DOS there. For 
U — 1 eV adopted in our calculation, the U c drops below 
U in narrow doping regions around the VH singularities, 
which will lead to magnetic long-range order. From the 
ferromagnetic correlation near the VH singularities re¬ 
vealed by shown in Figs. gc), ! e ), andgf), we 
conclude that long-range itinerant FM will emerge for 
U c < U in these narrow doping regions. 



X 

FIG. 4. The doping dependence of the magnetic critical inter¬ 
action strength U c - The horizontal solid line indicates U — 1 
eV, and the vertical dashed lines indicate the VH singularities. 


IV. TRIPLET p + ip' SC 


Away from the above introduced narrow doping region 
for itinerant FM, the interaction strength U is smaller 
than the critical value U c . Then through exchanging 
short-range spin or charge fluctuations between a Cooper 
pair, exotic superconducting states will emerge in the sys¬ 
tem as shown in Fig. [5] 

More specifically, we consider the scattering of a 
Cooper pair from the state (fe', — fe') in the /3-th (/3 = 
1, • • • ,4) band to the state (fe, —fe) in the <a-th (a = 
1, • • • ,4) band via exchanging spin or charge fluctuations. 
This scattering process leads to the following effective in¬ 
teraction vertex [ 34] : 

V aP (k, fe') = Re Y TP st ( fe > *')£* W 

pqst 

(io) 
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FIG. 5. The doping dependences of the pairing eigenval¬ 
ues A for all possible pairing symmetries. The vertical bold 
grey lines indicate the doping regions where the itinerant 
FM occurs. Inset: the typical split between the helical 
{px+ip y )n> ( Px-ip y )u pairing and the chiral {jpx±ip y )(n+m 
pairing caused by the weak SOC term with A so — 10 meV. 
The split in other doping regions where the p-\-ip' SC occurs 
is similar to the one shown in the inset. 


Here, for the singlet channel, we have 


r p s q t(k,k') = (U)%+\ [3(£/)(x (s) - x (c) m] Pt (k - k’) 

4 L J qs 

+ \ \3(U)(x {s) -X {c) m] PS (k + k 1 ), 

4 L \qt 

( 11 ) 

and for the triplet channel, we have 
C(fc, k’) = [(C/)(x (s) + X {C) )(U) f (k - k') 

4 L J qs 


(U)(x {s) + X {C) )(U)} PS (k + k’). (12) 
J qt 


From the effective interaction vertex ([Io| ), we obtained 
the following linearized gap equation |33j near the super¬ 
conducting critical temperature T c : 


1 



V a ^{k,k') 

Vp(k') 


A p{k') 


AA a (k). (13) 


Here the integration is along various FS patches labelled 
by a or /3, v^(k f ) is the Fermi velocity and fey is the com¬ 
ponent of k' along the FS. Solving this gap equation as an 
eigenvalue problem, one obtains each pairing eigenvalue 
A and the corresponding normalized eigenvector A a (k) 
as the relative pairing gap function. The leading pairing 
symmetry is determined by the A a (fe) corresponding to 
the largest A. The critical temperature T c is determined 
by A through T c = cutoff energy • e -1 / A , where the cutoff 
energy scales with the low energy bandwidth. 


According to its D^d point group, we study the pos¬ 
sible pairing symmetries of the system including s, p, d, 
and /-wave ones. The doping dependences of the largest 
pairing eigenvalues A for these pairing symmetries are 
shown in Fig. [5j At the low doping region, the doubly 
degenerate d x 2_ y 2 and d xy singlet pairings serve as the 
leading pairing symmetries, consistent with our previous 
results for the undoped case [22 . The gap function of the 
d x 2_ y 2 and d xy symmetries are symmetric and antisym¬ 
metric about the x-axis and y -axis respectively as shown 
in Figs. [6ja) and 0b). Below T c , these two degenerate 
pairing states will further mix to form the fully-gapped 
d x 2_ y 2 ± id xy (abbreviated as d + id') pairings to lower 
the energy. Physically, the d + id' singlet pairing is medi¬ 
ated by the antiferromagnetic spin fluctuations suggested 
by Figs. [3j(a), ^b), ^d), and^g). More importantly, 
in the vicinity of the narrow doping region for itinerant 
FM around the VH singularities, our RPA results iden¬ 
tify the doubly degenerate p x and p y triplet pairings as 
the leading pairing symmetries. The gap function of the 
Px{Py) symmetry is symmetric about the re(y)-axis and 
antisymmetric about the y(x )~axis respectively, with gap 
nodes on the ^(x)-axis, as shown in Fig. [6^c)([6][d)). The 
emergence of the triplet p x and p y pairings near the VH 
singularities is the physical consequence of the strong fer¬ 
romagnetic spin fluctuations there, as revealed by Figs. 
§c),§e), and^f). 


(a) | k y (b) f k y 
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FIG. 6. Distributions of the gap functions on the FS: (a) 
d x 2_ y 2 and (b) d xy symmetries for doping x = 0.12, as well 
as (c) p x and (d) p y symmetries for doping x — 0.24. 


Since the two p-wave pairing symmetries are degener¬ 
ate, they will probably mix to lower the energy below 
the critical temperature T c . To determine this mixture, 
we set A^ = Kip%(k) + (K 2 + iKs)p y (k ), where p x (k) 
and Py(k) denote the normalized gap functions of corre¬ 
sponding symmetries. Then the mixing coefficients K\, 
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K 2 , and K 3 are determined by the minimization of the V. DISCUSSION AND CONCLUSION 

total mean-field energy 


eg - A* _ 

v / (eg - M ) 2 + |A“| 2 _ 

+ 3^ E v “ l W) 

kk'a/3 

( A fc)* a£, 

Here the chemical potential p is determined by the con¬ 
straint of the average electron number in the supercon¬ 
ducting state. Our energy minimization gives K\ = ±K 3 
and K 2 = 0, which leads to the fully-gapped p x db ip y 
(abbreviated as p + ip') SC. This mixture of the two p- 
wave pairings satisfies the requirement that the gap nodes 
should avoid the FS to lower the energy. Note that there 
can be three different components of this triplet pairing 
with different S z quantum numbers of the Cooper pair, 
which are and (tl + it)- In the absence of SOC, 

the three spin components are degenerate. 

To lift up the degeneracy among the three different 
spin components of the triplet p + ip' pairing, we add the 
following Kane-Mele SOC term to Hamiltonian 0 : 

Hkm = ^A so ^ (15) 

<«» 

Here Vij = (2/v / 3)(hi x 62)2 = ±1 with b\ and 62 be¬ 
ing unit vectors along the two bonds that connect next- 
nearest-neighbors i and j on the same layer. Such a SOC 
term lifts up the degeneracy between the S z = 0 com¬ 
ponent and the S z = ±1 components. Our RPA calcu¬ 
lations (see the Appendix for the details) yield that the 
equal-spin helical (p x + ip y )(p x - ip y )u pairing wins 
over the chiral (p x =t ip y ) (fi+tf) pairing by a small split 
proportional to A so, as shown in the insets of Fig. [5] 
for A so = 10 meV. Such a helical triplet pairing is TR 
invariant weak topological SC. 

It is interesting to note that the pairing eigenvalue A 
of the p + ip' SC diverges in the doping region near the 
phase boundaries between FM and SC, due to the di¬ 
vergently strong ferromagnetic spin fluctuations in that 
region. Although the divergence of A is an artifact in 
the RPA caused by ignorance of the renormalization of 
the single-particle Green’s function, it is possible that 
the strong FM fluctuations in the critical regions con¬ 
siderably push up the superconducting critical tempera¬ 
ture T c , which might be experimentally accessible. Tak¬ 
ing into account that the doping level is hard to control 
in practice, we can instead apply tunable strain to the 
system to change the hopping parameters and the band 
structure [22]. As a result, the VH doping levels and the 
phase boundaries between FM and SC will shift so that 
a given doping level can access the phase boundaries to 
produce the high-temperature triplet p + ip' SC. 


The unconventional triplet p-\-ip' SC proposed here can 
be detected by various experiments. First of all, as an 
unconventional superconducting state with the phase of 
its pairing gap function changing on the FS, the p + ip' 
pairing state should show no Hebel-Slichter peak upon 
the superconducting phase transition in the NMR relax¬ 
ation rate 1/T\T [37] . Secondly, in this triplet pairing 
state, the Knight shift should not obviously change be¬ 
low the T c [38j. To further identify the phase structure 
of this pairing experimentally, we can fabricate a slice 
of BLS into a hexagon, and use a dc SQUID to detect 
the relative phase among different directions in the sys¬ 
tem [39] . In particular, determined by the p + ip'-wave 
symmetry, the phase difference between the opposite (ad¬ 
jacent) edges of the hexagon should be 1 r (tt/ 3). 

Although the p-wave SC is unconventional, the mix¬ 
ing of the p x and p y pairings into the complex p + ip' 
one leads to a fully-gapped superconducting state, which 
looks similar to the conventional s-wave one in many as¬ 
pects. For example, near zero temperature, the specific 
heat, the penetration depth, and the NMR relaxation 
rate in both fully-gapped pairing states decay exponen¬ 
tially with temperature. What’s more, the STM spectra 
of both fully-gapped superconducting states should ex¬ 
hibit U-shaped dl/dV — V curves. However, all these 
expected experimental results can be changed by a uni¬ 
axial strain applied on the system. More specifically, the 
p + ip' mixing proposed here is based on the degener¬ 
acy between the p x and p y pairing states, and the de¬ 
generacy itself originates from the D 3 d point group of 
BLS [ 22 ] . Thus, by applying a uniaxial strain to break 
the D 3 d symmetry of the system, we can eliminate the 
p + ip' mixing, and leave a single real p-wave pairing as 
the leading instability. Such a p-wave pairing can be the 
p x or p y one, which is determined by the axis of the ap¬ 
plied strain. Because the resulting p x or p y pairing has 
gap nodes on the FS, the above-mentioned exponential 
temperature dependence of the experimental quantities 
in BLS will be replaced by power-law ones. Meanwhile, 
the U-shaped STM spectrum of BLS will be replaced by 
a V-shaped one. 

In conclusion, we have systematically studied the pos¬ 
sible electronic instabilities of doped BLS. The results of 
our RPA calculations predict that the system is an itiner¬ 
ant ferromagnet in the narrow doping regions around the 
VH singularities, and a triplet p-\-ip' superconductor with 
a possible high T c in the vicinity of these regions. With 
an extra weak Kane-Mele SOC, we further single out the 
equal-spin helical p-\-ip' pairing state as the leading one. 
This intriguing triplet superconducting state has TR in¬ 
variant weak topological property, and can harbor the 
Majorana zero-mode at its boundary [3l I40H42] , which is 
useful in the topological quantum computation. 
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VI. APPENDIX: RPA WITH THE KANE-MELE 

SOC 


The Kane-Mele SOC term breaks the SU(2 ) spin- 
rotation symmetry, but keeps the C7(l) spin-rotation 
symmetry around the 5*-axis, leaving the ^-component 
of the total spin to be a good quantum number. There¬ 
fore, we define the following susceptibility tensors, 




k'a/3 


e < 22 ) 


A, , ( ~ ul ' I _ I 

lLJ n + £ k , 


xi t mpq (k,iu n ) = F E + k) 


k'a/3 

+ (23) 
lU) n + £ k r - £ k 7 +k 

x£ mpq (k,iu n ) = ^ E $(*')&/W&W + k) 

k' a/3 

&(>*'+< 24 ) 

lu} n + £ k > ~ £fc/+fc 

X { T )pq (k,iu n ) = ^(k'X^k^k' + k) 

k'a/3 

^(fc-+fc) " r(eg +»y n ^ ) . ( 25 ) 

iuj n + e: fe , 


In the RPA, we have 


xif P9 (fe,q 


(2 )pq 

Xsi 


(fc,r) 


xi? )p9 (fe,q 


(4)pqr 

Xst 


(fc,r) 


xif P9 (fe,q 


( 6 )pg 

Xlt 


( k,r ) 


v E ( TrC pt( fc i’ r ) c 9t( fe i + fc > r ) 

^ 1^2 

4 t (fc 2 + fe,0)ct t (fe 2 ,0)) , (16) 

v E ( T T c pt( fe i> r ) c 9t( fe i + fc > r ) 

^ 1^2 

c^(fe 2 + fc,0)c u (fc 2 ,0)), (17) 

jf E ( T r4l( fcl ’ r ) c 9l( fel + fc ’ r ) 

^ 1^2 

4 t (fc 2 + fc,0)c tt (fc 2 ,0)), (18) 

v E ( 7 v c pi( fei ’ r ) c 9i( fc i +fe ’ r ) 

kxk 2 

c h(k2 + k,0)cti(k 2 ,0)^ , (19) 

V E ( TrC pt( fcl ’ r ) c 9l( fel + fc ’ T ) 

fei k 2 

c^(fe 2 + fc,0)ct. t (fe 2 ,0)), (20) 

V E ( T ^ c pl( fe i’ r ) c 9t( fe i + fc , T ) 

fei k 2 


cl t (k 2 + fe, 0 )ct^(fe 2 , 0 )^ . 


X 


(5) 


X 


( 6 ) 


/ / x (i)(o) (E/q Vx (i)(o q 

h 4)(o y) / ) v ° ) 

, (26) 

( i x (i)(o) (^v7 0 i 
U (4)(o) (^) 1 ) U (4)(0 v 

, (27) 

J-X (5)(0) (t/) 

■■ x ( 5)(°) ; 

(28) 

J_ x (6)(0) (C/) ‘ 

_1 x (6)(0) , 

(29) 


where ([/) is the same as that in Eq. ©• 

With the above expressions of we consider the 

scattering of a Cooper pair from the state (fe', —k') in 
the /3-th ((3 = 1, • • • , 4) band to the state (fe, — k) in the 
<a-th (a = 1, • • • ,4) band. This scattering process leads 
to the following effective interaction vertices: 


qf (fe, fe') = Re E r ?u (fe, fe %7 (fe) 

pqst 

C(-fe)^(-fe , )ef t (fe'), (30) 

qf (fe, fe') = Re E rP tf (fe, fe')Ct* (fe) 

pqst 

ey(-fe)^t(-fe , )^t(fe')- 


( 21 ) 


(31) 










Here 


ing channels, 


pt 


?%l{k,k')=(u)$- (u)x {3 \u) (fe-fc') 




. . ips 

0 U)X {6 \U ) (fc + fe'), (32) 

J qt 


1 

2 L 

1 

2 L 


pt 


(U)x {4 \U) (fc-fe') 


J <?S 
ps 


(U)x {A \U) (fc + fc'). (33) 

J qt 


V“ e %(k,k') =- [vtf(k,k!) + vtf(k,-k') 

*£? 0) (*,*') [nf ( fc ’ fc/ ) - Pf ( fc ’- fe ') 

V^ ±l) {k,k') =V^(k,k'), 


(34) 

(35) 

(36) 


where the index e is for the even parity pairing, and o 
for the odd one. From these symmetrized effective inter¬ 
action vertices, we obtained the following linearized gap 
equation near the superconducting critical temperature 
T c : 


1 



^4)( fc ’ fc ') 

Vp(k') 


A p{k') 


XA a (k), 

(37) 


The inversion symmetry, together with the U( 1) spin- 
rotation symmetry of our system, enable us to sym¬ 
metrize the effective interaction vertices into the follow- 


which replaces Eq. (13) to determine the T c and the 
leading pairing symmetry of the system in the presence 
of the Kane-Mele SOC. 
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